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Abstract 

We prove an index theorem for the existence of Majorana zero modes in a semiconducting 
thin film with a sizable spin-orbit coupling when it is adjacent to an s-wave supercon- 
ductor. The theorem, which is analogous to the Jackiw-Rebbi index theorem for the zero 
modes in mass domain walls in one-dimensional Dirac theory, applies to vortices with odd 
flux quantum in a semiconducting film in which s-wave superconductivity and a Zeeman 
splitting are induced by proximity effect. The momentum-space construction of the zero- 
mode solution presented here is complementary to the approximate real-space solution of 
the Bogoliubov-de Gennes equations at a vortex core [J. D. Sau et al, arXiv:0907.2239], 
proving the existence of non-degenerate zero-energy Majorana excitations and the resultant 
non-Abelian topological order in the semiconductor heterostructure. With increasing mag- 
nitude of the proximity-induced pairing potential, the non-Abelian superconducting state 
makes a topological quantum phase transition to an ordinary s-wave superconducting state 
which no topological order. 



1 Introduction 



Particle statistics is a genuinely quantum mechanical concept which has no clas- 
sical analog. In spatial dimension three and above, the wave function of a many- 
body quantum state of identical particles remains either unchanged (bosons) or 
undergoes a change of sign (fermions) under a pairwise interchange of the parti- 
cle coordinates. That there are only two possibilities is a consequence of the fact 
that there are only two irreducible representations of the permutation group for N 
particles. However, in (2 + 1) dimension, where permutation and exchange are not 
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necessarily equivalent, the quantum statistics of particles can be remarkably differ- 
ent from the ordinary statistics of bosons and fermions HI 121 . In this case, under a 
pairwise interchange of the particle coordinates, it is possible for the many body 
wave function to receive an arbitrary phase factor e %e where 9 is an angle interme- 
diate between (boson) and n (fermion). The particles obeying statistics given by 
the angle 6 are called anyons [3J. Even if the anyon statistics is remarkably differ- 
ent from that of bosons and fermions, since the simple phase factor is only a one- 
dimensional representation of the braid group in 2D, the statistics is still Abelian. A 
more exotic possibility, one which has remarkable prospect for fault-tolerant topo- 
logical quantum computation (TQC) [4], arises when the quantum ground state of 
the many-particle system is degenerate. In situations where the many body ground 
state wave function is a linear combination of states from this degenerate ground 
state subspace, a pairwise exchange of the particle coordinates can unitarily rotate 
the wave-function in the subspace. In this case, the exchange statistics is given by 
a multi-dimensional unitary matrix representation of the 2D braid group, and, thus, 
the statistics is non- Abelian. It has been proposed [5] that such a system, where the 
ground state degeneracy is protected by a gap from local perturbations, can be used 
as a fault- tolerant platform for TQC. 

Recently, some exotic ordered states in condensed matter systems, such as the Pfaf- 
fian states in fractional quantum Hall (FQH) systems H6I7I8I9I and chiral p-wave 
superconductors/superfluids H10I11I12I13L as well as the surface state of a topo- 
logical insulator (TI) in which s-wave superconductivity is induced by proximity 
effect (1141151 1611 , have been proposed as systems which support quasiparticles with 
non-Abelian statistics, and, therefore, can potentially be used as TQC platforms. 
The common thread between these systems is that they all allow quasiparticle ex- 
citations which involve no energy cost (when the mutual separation among the 
excitations is large). The second quantized operators, 7,, corresponding to these 
zero energy excitations are self-hermitian, 7] = 7^ which is in striking contrast to 
ordinary fermionic (or bosonic) operators for which q ^ q. However, since 7* and 
jj, which are called Majorana fermion operators, anticommute when i ^ j, they 
retain some properties of ordinary fermion operators as well. It is the self-hermitian 
property of the Majorana operators which lies at the heart of the ground state de- 
generacy and the resulting non-Abelian statistics 1171 1011 of quasiparticle excitations 
in these systems. 

It has been shown recently [17J that even a regular semiconducting film with a 
sizable spin-orbit coupling, such as InGaAs thin films, can host, under suitable 
conditions, Majorana fermion excitations localized near defects. By an analysis of 
the real- space Bogoliubov-de Gennes (BdG) equations for a vortex in the semicon- 
ductor, in which s-wave superconductivity and a Zeeman splitting are proximity in- 
duced (Fig. (la)), it has been shown that the lowest energy quasiparticle excitation 
is a zero-energy Majorana fermion mode. The Majorana mode is separated by a fi- 
nite energy gap (so-called mini-gap) from the other conventional fermionic excited 
states in the vortex core. Thus, for a collection of well-separated vortices, the re- 
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suiting degenerate ground state subspace is protected from the environment by the 
mini-gap, enabling the potential use of the semiconductor hetero structure in Fig. 
(la) in TQC. Since the basic effects behind the emergence of the Majorana fermion 
excitations - spin-orbit coupling, s-wave superconductivity, and Zeeman splitting 
- are physically well-understood and experimentally known to occur in many solid 
state materials, the proposed semiconductor hetero structure ifTTll is possibly one of 
the simplest systems to realize, which supports non-Abelian topological order. 

In this paper, we show that the existence of the Majorana fermion zero modes at the 
cores of the vortices in the semiconductor hetero structure in Fig. (la) is due to an 
index theorem. The theorem is analogous to the one proven by Jackiw and Rebbi 
H18I19I for the existence of zero energy eigenstates at mass domain walls in a ID 
system of Dirac fermions, which was later applied ll20ll to prove the existence of 
zero-energy excitations at vortices in a spinless 2D chiral p-wave superconductor. 
In analogy with the chiral p-wave superconductor, we find that while an odd flux- 
quantum vortex in the semiconductor thin film traps a unique zero energy eigen- 
state, an even flux-quantum vortex does not. Furthermore, since the components of 
the fermion field in our analog of the ID Dirac theory are pairwise related by parity 
reversal and hermitian conjugation, we explicitly show that the second quantized 
operator corresponding to the zero-energy eigensolution is self-hermitian, 7^ = 7. 
This is in contrast to the corresponding problem solved by Jackiw and Rebbi for the 
conventional Dirac theory, where the two components of the fermion field are inde- 
pendent, and, therefore, the zero-energy excitations localized in the mass domain 
walls are conventional fermion excitations. The momentum space construction of 
the zero-energy solution presented here is complementary to the approximate real 
space solution |[T7l of the four-component BdG equations at a vortex in the semi- 
conductor in which superconductivity has been induced by the proximity effect. 



2 Hamiltonian 



The single-particle effective Hamiltonian H for the conduction band of a spin- 
orbit-coupled semiconducting thin film in contact with a magnetic insulator (which 
induces the Zeeman splitting) is given by (h — 1) H — H K + H S o + H z , where, 



k 2 

k,/3 z " t 

H so =-a^ [(k x a).z]p r/ cl tf3 c ktl + H.c. 

H Z =V Z J2 (^)/3,7 C k,/3 C k, 7 (!) 
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Fig. 1. (a): The proposed heterostructure of a semiconductor (SM) sandwiched between an 
s-wave superconductor (SC) and a magnetic insulator (MI). In this geometry, the semicon- 
ducting film can support non-Abelian topological order, (b): Single-particle band-structure 
in the semiconducting film with and without the Zeeman splitting. Without the Zeeman 
splitting, the two spin-orbit shifted bands touch at k = (red lines). Then, for any value of 
the chemical potential, the system has two Fermi surfaces. With a finite Zeeman splitting, 
the bands have an energy gap near k = (blue lines). If the chemical potential lies in 
the gap, the system just has one Fermi surface (indicated by the dotted circle), a situation 
conducive to the emergence of non-Abelian order. 

Here, m*, V z and [i are the conduction-band effective mass of an electron, effective 
Zeeman coupling induced by proximity to a magnetic insulator, and chemical po- 
tential, respectively. The Zeeman coupling can also be intrinsic if the film is made 
of a magnetic semiconductor (e.g., GaMnAs). We assume that the Zeeman coupling 
in the semiconductor is the dominant effect of the proximity to the magnetic insu- 
lator and the direct magnetic field induced by the magnetic insulator is neglected. 
The coefficient a describes the strength of the Rashba spin-orbit coupling in the 
semiconductor, /3, 7 are the spin indices, and a x , a y , a z are the Pauli matrices. The 
energy spectrum of the spin-orbit-coupled semiconductor with V z = has two 
bands crossing the Fermi level. This situation should be contrasted with the metal- 
lic surface state of a strong TI where an odd number of bands cross the Fermi level. 
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However, by placing the semiconductor in the proximity to a magnetic insulator, 
it is possible to open a gap in the spectrum, see Fig. (lb). Then, for |/i| < \V Z \, a 
single band crosses the Fermi level. This situation is then analogous to the surface 
of a strong TI in that an odd number of bands cross the Fermi level lfT4l . and, 
therefore, is suggestive of supporting non-Abelian topological order if s-wave su- 
perconductivity can be induced in the film. This way, by replacing the surface of 
a 3D strong topological insulator by a regular semiconducting thin-film, the scope 
of solid state systems which can be designed to support non-Abelian topological 
order can be vastly expanded IfTTll . 

In the presence of an adjacent s-wave superconductor, which induces an s-wave 
pairing potential in the semiconductor by the proximity effect, the full Hamiltonian 
in the bulk semiconductor becomes, Hb = H + H p , where, 

iJ p = -A ^4/ kii + H.c. (2) 

k 

We show below that the heterostructure in Fig. (la) has non-Abelian topological 
order by proving a theorem for the existence of a Majorana zero energy mode at 
the core of a vortex with an odd number of flux-quantum in the superconductor. Our 
momentum space construction of the zero mode solution in the form of an index 
theorem is complementary to the analysis of the real space BdG equations recently 
carried out IfTTl for the same system. 

In order to describe the spatial dependence of the superconducting order parameter 
in the presence of a vortex with winding number 1, we first write GUI the pairing 
part of the Hamiltonian in real space, 

H v = -Jd 2 Rj d 2 re ie ^h{R)g(r)(^ +r ^_ rA + H.c. (3) 

Here, R and r are the center-of-mass and the relative coordinates of a Cooper pair, 
respectively. h(R) and 6>r, are the amplitude and the phase of the superconducting 
order parameter, and, in momentum space, g(k) is the s-wave order parameter A . 

In the vortex core, h(R) ~ (1 — e~ T) with £ the coherence length, which describes 
the suppression of the condensate amplitude at the core of a vortex. In order to 
prove the existence of a non-degenerate zero energy state at the vortex core, we 
have to map the 2D Hamiltonian H T = H + H v on the Jackiw-Rebbi problem of 
ID Dirac theory with a mass domain wall. Before we do this, let us briefly discuss 
the derivation of the zero-mode in the Jackiw-Rebbi problem [fT8l itself. 
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3 Fermion zero mode in ID Dirac theory 



Let us briefly review the derivation of the Jackiw-Rebbi zero mode for the ID 
problem described by the Dirac Hamiltonian, 



Hn = dx 



iv Fif) cr z d x ip + ■m(x)i/j^ a x ijj 



(4) 



where ift(x) = (fl(x), flfa)) with f\,2{x) two independent fermion fields. Here, 
vf is the Fermi velocity of the fermions and m(x) is a spatially- varying mass for 
the fermion fields. To uncover the form of the real space wave functions of the 
quasiparticle excitations, we first write the second-quantized quasiparticle operator, 

qi = Jdx [MtffRx) + Ux)fl{x)l (5) 
which is assumed to satisfy the equation, 

[H,q*] = eq\ (6) 

This implies the following real-space Dirac equation for the two component wave 
function (j) T (x) = (4>i(x), (f) 2 (x)): 

— ivF(y z d x <^{x) + a x m(x)4>(x) = e<p(x). (7) 

As a first step towards solving this equation, we note that, because a y anticommutes 
with a x and a z , if <f>(x) is an eigenfunction of Eq. © with eigenvalue e, a y (j)(x) 
is also an eigenfunction of Eq. © with eigenvalue — e. As a result, if there is a 
non-degenerate solution of Eq. © with the eigenvalue e = 0, it can be made a 
simultaneous eigenstate of a y since a y 4>(x) oc (j)tx). Let 4>q{x) denote such a zero- 
energy solution and a y 4>o(x) = X<po(x). Now, setting e = and left-multiplying 
Eq. © by ia z we obtain, 

d x <fi (x) = —m(x)<j>o(x), (8) 
v F 

which implies 

Mx) = e^ I ° m(y)dy M0)- (9) 
It can be clearly seen that for m(x) = ±sign(x) \m(x) |, Eq. © is normalizable for 
A = =Fl. Therefore, for each sign change of the mass term m(x) in the ID system, 
there is a single normalizable zero energy solution. Such a solution indicates the 
existence of a zero-energy excitation with a quasiparticle operator (for A = 1), 

t = C j dxe^^ m{v)dy [fl{x)+ift{x)l (10) 



Q 

where C is a normalization constant. Since fi(x) and f%(x) are two independent 
fermion fields, above denotes an ordinary fermion operator following fermion 
anticommutation relations. 
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Below we map the 2D Hamiltonian Ht = H + Hy of the semiconducting thin 
film with proximity-induced superconductivity on an effective ID problem which 
resembles Eq. ©, with important differences which render the quasiparticle mode 
corresponding to the zero-energy solution a Majorana fermion mode. 



4 Map of the Hamiltonian of the 2D semiconductor on an effective ID theory 



In order to map the Hamiltonian of the semiconductor, H T = H + H v , with 
proximity-induced superconductivity and a vortex at the center, on an effective 
one-dimensional problem, we will use the rotational symmetry of H T to decou- 
ple it in the various angular momentum channels EOl . We note that, because of 
the existence of the spin-orbit coupling term H so in H , the Hamiltonian H T is 
invariant only under the simultaneous rotation of the system in the real and spin 
spaces. Therefore, for angular momentum decoupling of H T , we change to a new 
representation of the fermion operators where they are expanded in the total an- 
gular momentum channels (as opposed to the orbital angular momentum channels 
as employed in the case of a chiral p-wave superconductor EOl ) indexed by the 
half-odd-integers mj, 



1 oo 

c . t = 1 V c u fe i{mj -^ 9k 

V t mj =- 00 
1 oo 

c *l=^TT £ c mjAl e^ + ^ d \ (11) 

V^TTfc mj= _ 00 

where (m j =)= 1/2) in the exponents on the right indicate the orbital angular momen- 
tum quantum numbers for the up and the down spins, respectively. Here, k = |k| 
is the magnitude of the vector k, and thus is a one-dimensional variable. The com- 
mutation relation {c k/3 , cj, } = 5 2 (k — p)<5g, 7 implies 

{Cmj,k,(3, 4ij,p,y} = $mj,nj5(k - p)5 0>T (12) 

Inserting these representations in Eq. O, we find for H K , 

Hk = pEE/*(^-^^W- d3) 

The Zeeman term, Hz, is also diagonal in the total angular momentum as well as 
the spin quantum numbers, 

V f 

Hz = J dk (4njM C ™J,kA ~ C Lj,k,l C mj,k,l)- ( 14 ) 



mj 
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The spin-orbit term, Hso, is diagonal in the total angular momentum quantum 
number, but not so in the spin quantum number: 

Hso = J kdkc !nj,k,1 C ™jAl + H.c. (15) 



Let us now focus on the vortex term H v in Eq. ©. Substituting 

c R±r,/3 _ Z7r C k,/3 e 
k 

in Eq.[3l we end up with two spatial integrals, g(k — p) = / <i 2 rg(r)e^ k ~ p ) r = A 
and 

J(k + p) = J d 2 Re i0R h(R)e iik+p) - R . (16) 

In order to evaluate I(k + p), we first note that, if in Eq. (fT6l) the vector (k + p) 
is rotated by an angle 9 in the momentum space, the scalar product in the exponent 
exp[i(k + p).R] can be made to remain invariant if simultaneously 9 R is rotated by 
9. Therefore, 

J( J R,(k + p)) = e ie /(k + p), (17) 

where Re is the operator that rotates (k + p) by an angle 9 in the momentum space. 
It follows that, I(k + p) = e iek +p/(|k + p|). To evaluate I(\k + p|) we choose 
(k + p) along the y-axis. Performing the #r integral which produces — 2ni J_i(|k + 
p\R), where J_i is the Bessel function of the first kind of order —1 |2T]|. and then 
performing the R integral which produces ij^p x O(l), we find, 

H V = -(2n)HA ^ ,^ P l 4, T ^,| + H - c - (18) 

I T P| 



k,p 



Finally, using angular momentum expansion of the fermion operators and noting 
that a function of |k + p| is periodic in (0 k — 9 P ) and hence can be Fourier expanded 
as 

' Y, u m(k,p)e m ^~ dp) (19) 



|k + p 
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with m an integer, we find, after the 9 k and 9 P integrals, 



(20) 

From Eqs. (1131 141 1 51201) . it is clear that the mj = 1/2 total angular momentum 
channel separates from the rest (in Hy it is obtained by taking m = —1 in the 
first term and m = in the second term on the right hand side of Eq. (|2Q|)). It is 
straightforward to check that no other total angular momentum channel is isolated 
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from the rest. One can also check that for a vortex with an even number of flux 
quanta, for which e l9R is replaced by e 2l9R in Eq. ©, no total angular momentum 
channel can be isolated from the rest of the Hamiltonian. The fact that a single total 
angular channel can be isolated from the rest is crucial for the following arguments 
demonstrating the existence of a non-degenerate zero energy solution. The zero 
energy eigen- solution of the Hamiltonian exists in this channel, and in cases where 
such an isolated channel does not exist, no non-degenerate zero energy solution 
exists at the vortex core. 

In the rrij = 1/2 channel, the total Hamiltonian can be written as H T \ = 

H K,m J= \ + H Z,m J= \ + H SO,m J= \ + H V,mj=± > where > 



= j^y^I dk ^ -ti c W\** 

iat r | 
H so,mj=\ = 7^31 J kdkc ±,k,f^M + H ' c - 

H v, mj =\ = ~i J dkdpA(k,p)c[ k ^c\ p l + }I.c., 

(21) 



where, in the last line, A(k,p) = A y/kp(ku m= -i(k, p) + pu (k,p)). From here 
onwards we will ignore the subscript 1/2 in the fermion operators c, c\ keeping in 
mind that all the fermion operators carry the same total angular momentum quan- 
tum number mj = 1/2. 



5 Demonstration of the zero energy solution 



To turn the Hamiltonian in Eq. (12TT) into a form resembling that in Eq. ©, we have 
to first diagonalize the free-electron part, (H Kmj= i +H Zmj= i +H SOmj= i), using 
a unitary transformation. Only then it is possible to linearize the band energies 
around the Fermi surfaces and the combined one-electron part in Eq. (|2~TT) may take 
the Dirac form. To do this, we define the unitary transformation from the cj n k 
basis to the energy eigenbasis, 
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\f-,k) 




(22) 



(23) 



where |afc| 2 + \b k \ 2 = 1. Since the fc operators are related to the original fermion 
operators by a unitary transformation, they still obey the fermionic anticommu- 
tation relations. Calling the corresponding eigenvalues (band energies) E± >k , the 
Hamiltonian in this basis in the channel rrij = | is written as, 



/elk 



(27T) 



-i / dkdpA(k,p)[alblfl k fl jP - b k a p f[ k f j _J 

- i J dkdpA(k,p)[a* k a p fl k f_ p - b k b* p f[ k flj + H.c. 



7 dMp[{A(A ;) p)46;-A(p ) A ; )a;^}4 )fc 4 jP ] 
| y dkdp[{-A(k,p)b k a p + A(p, k)b p a k } f ] _ k f ] _^ p ] 

1 I dkdp[{A(k,p)ala p + A(p,k)b p bl}fljlj+E.c. 



(24) 



Here, in the second expression, we have written the coefficients of the pairing terms 
in a form so as to highlight the fact that the coefficients of the intra-band pairing 
terms must be antisymmetric under the interchange of k and p, while the coeffi- 
cients of the inter-band pairing terms have no such symmetry. The antisymmetry of 
the coefficients of the intra-band pairing terms, which will be crucial for the demon- 
stration of the zero-energy solution below, follow from the fermion anticommuta- 
tion relations coupled with the fact that the fermion operators involved are from the 
same band with the same mj quantum numbers. Using simplified notations for the 
coefficients of the pairing terms we rewrite Eq. (|24l) as, 



/fl rC 

2J dkdp[{A ++ (k,p)fljl tP + A__(fc,p)/V-.p 
+ A + 4k,p)fljlj}+R.c. 



(25) 



Let us now assume that the Fermi level is in the lower band /_, which has a Fermi 
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momentum kp, and that the upper band / + has a set of unoccupied states at energy 
E +>k — n separated from the Fermi level by an energy gap E + = \V Z \ — fx (without 
loss of generality we take fi to be positive). This condition is satisfied only as long 
as I/j-I < \V Z \. In what follows we will implicitly assume that this condition is 
experimentally satisfied. By limiting the momenta in the lower band to those close 
to the Fermi momentum and in the upper band to those close to the band minimum 
at k = 0, the Hamiltonian can be written as, 



Ht = J ^^/U+tS-,^,-/!)/*^,/- 



,k F +q\ 



- l -J dqdq'[{A ++ (q, q')fljl q , + A__(A* + q, k F + ffifl^ft^ 
+ A+_(g, k F + q')fljl kF+q ,}] + H.c. (26) 



Here, the lower limits on the momenta q, q' in the integrals are or —kp when the 
corresponding momentum corresponds to the + or the — bands, respectively (the 
upper limits on the integrals are taken to be oo). We keep the lowest order depen- 
dencies of the functions A on q and q' which are consistent with the symmetries of 
these functions from the fermion anticommutation relations, namely that A ++ and 
A are antisymmetric in (q — q'). Since A + _(q, q') has no such symmetry require- 
ment, it is generically non-zero for q, q' — > 0, and we assume that it is a non-zero 
constant A + _ in this limit. This yields, 



Ht = J j2^[E + fl jq f +jq + (E^ ikF+q - v)f! >kF +qf-,k F +Q 



J dqdq'[{A ++ (q - q')fl >q fl, qt + A-(q - q')fl, kF+q fl 



k F +q' 



+ A + _/i, 9 /i fcF+9 ,}]+H.c. 



(27) 



Let us now define the Fourier transforms to real space, 



fl{x)= r dqe^fl kF+q 

J —kp 

fl(x)= r°dqe^fl A . (28) 

J 



Using Eq. (|28l) , we fourier transform H T to real space, 
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H T = r dxdx'[E + 5(x - x')fl(x)f + (x') + E^(x - x')fl(x)f-(x')} 



oo 



-~j dx[{A ++ (x)fl(x)fl(-x) + A__(x)/1 (x)fl(-x) 

+ A + _8(x)fl(x)fl(-x)}}+E.c, (29) 



where the functions A ++ (x) and A (x) are odd in x, and E-(x—x') = ivd x 5(x — 

x') with v the fermi velocity in the lower band. The Hamiltonian in Eq. (|29l) is still 
not in the same form as in Eq. ©, so the Jackiw-Rebbi construction does not yet 
apply. However, as we show below, the BdG equations for a putative zero-energy 
excitation in the vortex core are of the same form as in Eq. © (with e = 0), 
guaranteeing the existence of a zero energy solution. 

Let us define the quasiparticle operator for the putative zero energy mode as, 



7 f = / dx(rj +tl (x)fl(x) + r}-x(x)fl(x) + rj +t2 (x)f+(-x) + T]^ 2 (x)f-(-x)). 

(30) 



The corresponding BdG differential equations can be derived by setting 

[#r,7 f ]=0. (31) 



To derive the commutators involved in Eq. (I3TI) . we need the following anticom- 
mutation relations, 



{A(x), /+(*')} = / dkdk'e< k *- k '^{fl >k , f +>k ,} 
Jo 

= / dke lk{x ~ x) =S+(x-x') 
Jo 



■ 5{x-x')/2 + iP{—^—;) (32) 



{/!(*), Ux')} = / dkdk'e^- k ' x '\fl k , 

/oo 
dk( ?k{x-x) = S _(x-x'), (33) 
-kp 



where the functions S + and SL are not exactly 5 functions. The existence of S + 
and S- necessitates the introduction of a new set of four functions, 



iv,n{x)=i S u {xi - x)r) v , n (xi)dxi, (34) 

J —oo 

where v = ± and n takes the values 1,2. 

In terms of the £ functions, the BdG equations take the simple form, 
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ivd x Z- A (x) - iA__ (x)£- 2 (x) -?A + _5(x)£ + , 2 (»/2 = 
-*A ++ (x)£ + , 2 (x) + *A + _<5(x)£_, 2 (x)/2 = 
+iA__(a:)^_,i(x) - 2A + _5(x)£ +i i(x)/2 = 

-£ + £ + ' 2 (z) + 4 A ++ (z)£+'i(z) +iA+_(J(a;)^i(x)/2 = (35) 

Eliminating the amplitudes in the upper band, £ +)2 (x), using 



tA + „5{x) I 
' 2(E 2 + + A 2 ++ (x)) ,a 



-iA. 



{x 



-E, 



£-,i(*) 



(36) 



we get the equations for the amplitudes in the lower band, 



1 A 2 5 2 {x)E + , , 4 



1 A^__(5 2 (x),E + 



UP 



1 A^_5 2 (x)A. 



UP 



z A fx) 



4 £ 2 +A| + (x) 
1A 2 + J 2 (i)A ++ (j;) 
4 M + A3..(x) 



= o 

= 
(37) 



It can now be easily seen that by a redefinition, 

_ _l r x d , A l- s2 ( x ' )E + 
M*) = e~ l4vJ ° £-,„(*), (38) 

the above BdG equations can be recast as, 



vd x ipi{x) + \A—(x) + 



lA 2 + _5 2 (x)A ++ (x) 
4 E\ + Al + {x) 



o , / n /a , s ^A 2 + _5 2 {x)A ++ {xY 
vd x ip 2 (x) + A__(x) + -- 1 



4 £ 2 +A 2 ++ (x) 



i[) 2 (x) 
ipi(x) 








(39) 



By defining the quantity in the parentheses in Eq. (l39l) as m(x), it can be easily 
checked that m(x) is an odd function, m(—x) = —m(x), since both A ++ (x) and 
A (x) are odd under x — > —x. Eq. (1391) can now be written in terms a two- 
component spinor wave function ( T (x) = (4>i(x), ^{x)), 





= -m(x)a x 


<^{x)\ 









(40) 



which, after multiplication by a z , is of a form similar to Eq. ©. Therefore, a unique 
zero energy eigen- solution of H T in the rrij = | channel is guaranteed. 
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Notice that, in the above construction for the zero-energy solution of the BdG equa- 
tions, no special form for the real-space profile of the superconducting gap near the 
vortex core has been assumed. In particular, we have not assumed the frequently- 
used step function profile for the gap [| 1 31 1711 near the core, without which the real 
space solution of the second-order, coupled, BdG equations for the present system 
is still an open problem. All we have needed here to prove the existence of the zero- 
energy solution, irrespective of the real space profile of the gap, is the antisymmetry 

of the functions A ++ (x) and A (x) under x — > —x, which is a consequence of the 

fermion anticommutation relations. Note that, to use the fermion anticommutation 
relations to determine the symmetry of these functions, it is crucial that a single to- 
tal angular momentum channel (rtij = 1/2) could be isolated from the rest. In the 
cases where this is not possible, e.g., a vortex with an even number of flux quanta, 
the Jackiw-Rebbi construction does not apply and a non-degenerate zero-energy 
solution is not expected at the vortex core. 

It is worth pointing out that the above construction for the zero-energy solution 
can be straightforwardly applied to describe the zero modes on the surface of a 3D 
strong TI, described by a Dirac-like Hamiltonian, in the presence of an s-wave su- 
perconducting vortex H14I22H . but only for a non-zero chemical potential. For zero 
chemical potential, our reduction of the set of 4 coupled BdG differential equations 
to a set of 2 (with the help of Eq. (1361 )) to bring them to the form in Eq. CZ]) in 
terms of 2 x 2 Pauli spin matrices no longer applies. Therefore, our methods above 
do not directly apply to the problem of zero modes in a TI in the presence of an 
s-wave superconductor at fi = 0. However, for this special case, there exists exact 
real space solution of the zero-mode eigenf unction H14I22H . 



6 Demonstration of the Majorana condition 

Now we show that the quasiparticle operator for the zero energy solution found 
above satisfies the Majorana condition 7* = 7. From Eq. (|30l ) it is clear that such a 
condition is satisfied provided 




—x) = 77- 2 (x) 

~X) = T] +j2 {x) 



(41) 



Using Eq 



(1341) and the property for the S± functions, 



S±{x\ — x) = S±(x — x\) 



(42) 



we see that the condition in Eq. (14TT) translates into 
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Now, from Eq. (|40l) it is clear that ((x) can be taken as an eigen-spinor of a x : 

a x ((x) = \((x) (44) 
Then, solving for Eq. (|4Q|) . we find that, 



A r x 

i)i(x) = exp( / m{y)dy) 

v Jo 

A r x 



a r x 

if) 2 {x) = sign(A) exp( / m(y)dy), (45) 

f Jo 

where, for m(x) = ±sign(x) \m(x) |, A = ±1, in which case the zero mode solution 
is normalizable. Here we take A = 1 for the purpose of illustration. Using Eq. (1381) , 
we get the solutions for the functions £_ i, £ 2 > 



1 r x A 2 _5 2 {x')E + A r x 

£_ A ( x ) = exp(z— / dx' J 2 A )exp( — / m(y)dy) = £_ 2 (a0 (46) 



Now, with the help of the phase factor in Eq. (1461) being even in x and the mass m\x) 
being odd in x, it can be easily checked that the first line in Eq. (|43l) is satisfied. 
Subsequently, with the help of Eq. (f36|) . and recalling that A ++ (x) is odd in x, 
the second line of Eq. (1431) is also satisfied. Therefore, the quasiparticle operator 
corresponding to the zero mode solution in Eq. (|46|) is a Majorana fermion operator. 



7 Topological protection and topological quantum phase transition 



The proof of the index theorem presented above relies on the rotational invari- 
ance of the non-interacting part of the Hamiltonian as well as of the vortex profile, 
together with the assumption that a single band (lower band, /_) dominates the 
wave-function of the Majorana bound state. The latter assumption is valid in the 
limit E + — \Vz\ — ^> A . However, as shown below, the question of the exis- 
tence of the non-degenerate Majorana mode in this system is robust to deviations 
from such assumptions, and the Majorana state can only be removed by a bulk 
phase transition in the system. Thus, a non-degenerate Majorana bound state is 
topologically protected. 

In order to see the topological protection of the Majorana mode, let us consider an 
isotropic BCS Hamiltonian Hi with a symmetric vortex and a small value of A 
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(basically, Hi is the same as Ht), for which the derivation presented above can be 
used to show the existence of a non-degenerate Majorana mode bound to a vortex 
with an odd number of flux quantum. In addition, let us consider a second BCS 
Hamiltonian H 2 which is possibly anisotropic {i.e., breaks rotational symmetry) 
near the vortex and has a strong pairing potential. Therefore, H 2 does not lend 
itself to a simple analytic treatment as given above. Given these two Hamiltonians 
H\ and H 2 , we can construct a family of BCS Hamiltonians parameterized by < 
A < 1 given by H(X) = XHi + (1 — \)H 2 . Since both Hi and H 2 approach the 
bulk Hamiltonian away from the vortex core, so does H(X). Therefore, as long 
as the superconductor has a non-zero bulk gap, E g , away from the vortex core, 
we expect all the eigenstates of H(X) with energy less than E g to be localized at 
the vortex core. Furthermore, since only a finite number of states can be localized 
near the vortex core, there should be a discrete set of eigenstates with energies 
satisfying \E^ X )+n)(X)\ < E g bound to the vortex core. Here, in analogy with a 
p-wave superconductor/superfluid ||23l , (s(X) + n) is an index for the states with 
E( s (\)+n) oc s(A) + n, where n is an integer and s(A) = 0, 1/2. Note that, if for 
some value of A, s(A) = 0, a single zero energy state is allowed. Conversely, if 
s(A) = 1/2, no zero mode is allowed. In general, the shift s(A) can be taken such 
that it represents an even number (including 0) zero modes by s(A) = 1/2 and an 
odd number of zero modes by s(A) = 0. The quasiparticle excitations of H(X) are 
given by the equation, 

d\ s{x)+n) (X)\ = E {s{x)+n) (X)d\ six)+n) (X). (47) 

where d\d are taken as the quasiparticle creation and annihilation operators. By 
considering the Hermitian conjugate of the above equation, it can be seen that the 
solutions of H (A) are paired such that d j _ (s(x)+n) (A) = d {s{x)+n) (X) and£_ (s(A)+n) (A) = 
— i?( s (A)+ n )(A). In the case where H(X) has a single zero energy solution, (s(A) = 
0), it is clear that the corresponding zero mode d is a Majorana mode (d\ = d ). 
Since H 1 in our problem is assumed to have a non-degenerate Majorana mode 
(as we have shown by analyzing H T ), s(X = 1) = 0. The stability of this non- 
degenerate Majorana zero mode to small perturbations in the BdG Hamiltonian 
follows from the fact that as long as the low energy states in the vicinity of the 
vortex core remain discrete in energy spacing, the eigenstates labeled by the index 
s(A) + n evolve continuously with A, and thus the shift s(A) cannot jump from it's 
value to it's other allowed value 1/2. Thus, s(0) = 0, which implies that H 2 has a 
non-degenerate Majorana mode as well, even though H 2 itself does not lend itself 
to a simple analysis. 

From the above argument it is clear that a non-degenerate Majorana bound state at a 
vortex remains protected as long as the low energy states in the vicinity of the vortex 
core remain discrete in number for all values of the parameter < A < 1. Thus the 
Majorana bound state is robust to all local deformations of the Hamiltonian in the 
vicinity of the vortex core since these cannot affect the energy gap away from the 
vortex. On the other hand, if changing the parameter A changes the bulk parameters 
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of the Hamiltonian, such as a, \x or V z , then it is possible for the energy gap to close 
away from the vortex core, causing the states bound to the vortex core to become 
delocalized. This results in a breakdown of the above argument for the topological 
stability of the Majorana bound states. Thus, the presence or absence of a non- 
degenerate Majorana bound state at a vortex core leads to a classification of the 
phase diagram of two-dimensional superconductors with spin-orbit coupling, such 
that the phases which support non-degenerate Majorana bound states are separated 
from the phases that do not by a quantum phase transition (QPT) where the bulk 
energy gap closes. 

The quantum phases on the two sides of the above QPT can be determined by 
studying a deformation parameter A that controls the proximity-induced pairing 
potential in the semiconductor, i.e., the pairing potential in H(X) is given by AA . 
The pairing potential increases from (A = 0) to a maximum value A (A = 
1). The proof presented above shows that for < \Vz\, for which there is a 
single band crossing the Fermi-level, and for AA <C E + = (\Vz\ — (jl), there 
is a non-degenerate Majorana mode bound to the core of a vortex. From the bulk 
Hamiltonian H(X), it is easy to check that the energy gap in the bulk vanishes at a 
critical value 1171241 . 



Thus our argument for the topological protection of the Majorana bound state im- 
plies that a Majorana bound state exists at the vortex core in H(X) as long as A < A c 
or A < A c . Moreover, this phase is separated from the phase with A > A c by a 
QPT at which the single particle energy gap vanishes. The spin-dependent terms 
in the Hamiltonian, such as the Zeeman splitting and the spin-orbit coupling, cease 
to play a critical role in the phase A > A c , and it is possible to reduce both these 
couplings to zero without crossing another phase transition (i.e., without E g going 
through again). Thus, the phase A > A c must be in the same class as a conven- 
tional s-wave superconductor without spin-orbit coupling or Zeeman splitting, and 
therefore cannot support a non-degenerate Majorana mode. Furthermore, the sym- 
metries of the superconducting order parameters are identical in these two phases. 
Therefore, these two phases, distinguished by the presence or absence of a Majo- 



rana mode, are separated by a topological QPT at A = JV 2 — ji 2 where the energy 
gap closes. 



8 Conclusion 

We prove a theorem for the existence of Majorana zero modes in a semiconduct- 
ing thin film with a sizable spin-orbit coupling, in which s-wave superconductivity 
and a Zeeman splitting are induced by proximity effect (Fig. (1)). Our momentum- 
space construction of the zero-mode solution in the form of the Jackiw-Rebbi in- 
dex theorem ||T8l , which is complementary to the approximate real- space solution 




(48) 




17 



of the BdG equations at a vortex core ifTTll . proves the existence of non-degenerate 
Majorana fermion excitations localized at the vortices in the semiconductor het- 
erostructure shown in Fig. (la). 

In our method, no special form (in particular, the frequently-used step-function 
form H13I17I0 for the real-space profile of the superconducting pairing potential 
near the vortex core is required. We use only the antisymmetry property of the 

intra-band pairing potentials, A ++ (x) and A (x) (Eq. (l29l )) under x — > —x, which 

is a consequence of the fermion anticommutation relations. Therefore, the non- 
degenerate Majorana mode found here is robust to local spatial deformations of 
the order parameter. For the construction in this paper to apply, it is crucial that a 
single total angular momentum channel (for a vortex with a single flux quantum: 
mj = 1/2) can be isolated from the rest of the Hamiltonian. In the cases where this 
is not possible, e.g., for a vortex with an even number of flux quanta, the theorem 
does not apply and a non-degenerate zero-energy solution is not expected at the 
vortex core. The methods of this paper can be straightforwardly applied to the zero 
modes on the surface of a TI with a non-zero chemical potential in the presence of a 
superconducting vortex IfLfll . but not to the case when the chemical potential on the 
TI surface is zero. However, in the latter special case of zero chemical potential, 
an exact solution of the BdG equations already exists H 141221 for the zero-mode 
eigenfunction for arbitrary real- space profile of the order parameter. 
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